function v = wtp_function(alpha,v_grid,c_with,c_without)
f = zeros(size(v_grid));

for i = 1:length(v_grid),
    f(i) = cara_function(v_grid(i),alpha,c_with,c_without);
end

[f_min index] = min(f); 
v = v_grid(index);	